The average shape of a fluctuation: universality in excursions of stochastic processes 
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We study the average shape of a fluctuation of a time series x(t), that is the average value 
(x(t) — x(0))t before x(t) first returns, at time T, to its initial value x(0). For large classes of 
stochastic processes we find that a scaling law of the form (x(t) — x(0))t = T a f(t/T) is obeyed. 
The scaling function f(s) is to a large extent independent of the details of the single increment 
distribution, while it encodes relevant statistical information on the presence and nature of temporal 
correlations in the process. We discuss the relevance of these results for Barkhausen noise in magnetic 
systems. 
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Many experimental investigations of a wide range of 
complex systems involve the measure of some scalar ob- 
servable over long time intervals, during which the sig- 
nal exhibits nontrivial fluctuations around some average 
value or avalanche-like bursts of activity separated by 
quiescent intervals A few examples are Barkhausen 
noise in magnetic materials |2(, solar flares in astro- 
physics [3, seismic activity in geophysics Q, or prices 
in financial markets The statistical features of such 
fluctuations reflect the properties of the dynamics that 
generates them, and their analysis is a key point for un- 
derstanding the system under investigation. Traditional 
tools for characterizing temporal series are correlation 
functions, distributions of durations of fluctuations and 
of their amplitudes. Here we focus on the average shape 
of a fluctuation, showing that it contains crucial pieces of 
information about the nature of the underlying process. 
Since this quantity can be easily extracted from experi- 
mental data, it is an useful tool in the statistical analysis 
of temporal series. 

The average shape of a fluctuation has been lately the 
subject of interest in the study of the Barkhausen ef- 
fect 0, 0, > where it has been used as a stringent test 
for the validity of theories against experiments. All cur- 
rent models, successfully used to describe most of the 
features of the phenomenon, are however unable to re- 
produce the form of fluctuations (pulses) measured ex- 
perimentally [8|. In particular, while all models predict 
a shape symmetric with respect to its middle point, em- 
pirical data yield a skewed form. This indicates that 
the understanding of the Barkhausen noise is still incom- 
plete. Moreover it shows that the average shape of a 
fluctuation is a sharper tool for discriminating universal- 
ity classes than critical exponents . Motivated by these 
recent results, we approach theoretically the problem of 
finding the average shape of a fluctuation for a generic 
signal, and of understanding which statistical features 
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FIG. 1: Schematic illustration of the definition of the average 
shape of a fluctuation. Above, the signal x(t) is shown and 
two successive positive fluctuations above x(0) are indicated. 
Below the two fluctuations are rescaled and compared to the 
average shape. 



of the process generating the signal are encoded in this 
curve. 

Let us define precisely what we intend by average fluc- 
tuation of a generic signal (Fig. QJ. If x(t) is the am- 
plitude of a signal, we define as duration T of the pulse 
the time needed for x(t) to return for the first time to 
x(0). Collecting positive pulses of the same duration we 
introduce the average amplitude of the signal after time 
t elapsed from the beginning of a fluctuation 

(x(t)-x(0)) T . 

Interpreting x(t) as the position of a one-dimensional 
walker at time t, (x(t) — x(0))t is the average trajectory 
of the walker before it returns to the starting position 
x(0), and is obtained by summing over all positive walks 
starting at time t = in x = x(0) and constrained to 
return for the first time to x(0) at time t = T. 

In this Letter we compute this average trajectory for 
stochastic processes of the form 



x(t + l) = x(t)+£(t) 



(1) 
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with x(0) =0 0- The increment is a random vari- 
able. Despite the simplicity of the processes considered, 
we uncover a rich phenomenology. We find that, in all 
cases considered, a scaling law is obeyed 



(x(t)) T =T a f(t/T), 



(2) 



in appropriate time regimes. The value of the exponent 
a results to be equal to the wandering exponent of the 
"free" process (i. e. not constrained to return to at 
time T, nor to be always positive). The scaling function 
f(s) turns out instead to be a sensitive probe of temporal 
correlations of the dynamical process, while it is to a large 
extent independent on the statistics of single increments. 
In particular, for uncorrelated noise the scaling function 
is always proportional to a semicircle. For correlated 
noise the tails are universal with respect to the single 
step distribution, and they depend only on the short- or 
long-ranged nature of the correlations. 

We consider processes of the form and focus our 
attention on the trajectory of a walk originating in x(0) = 
at time and returning to x(0) for the first time at time 
T, known in the mathematical literature as excursion. 
We first study the case of uncorrelated random walks 
and Levy flights and then move to the analysis of the 
effect of correlations. 

i) Gaussian walks: The simplest process of this kind is 
the uncorrelated unbiased random walk, where (£(£)} = 
and = $t,t'- By using the time-reversal sym- 

metry the average trajectory can be written as 



r o °° dx x F{x,t)F(x,T -t) 
f^dx F(x,t)F(x,T -t) 



(3) 



where F(x, t) is the probability that a walker starting in 
x = at time reaches x at time t without ever touching 
the axis x = 0. 

If we consider a normal distribution for single steps 
F(x,t) can be determined via the image method |10| . 
The result for (x(t)) T is of the form (0) with a = 1/2 
and scaling function proportional to a semicircle 



f u ( S ) = x l-^r^7). 



(4) 



The exponent a, being equal to the value characteriz- 
ing a free walk, indicates that the constraint of returning 
at time T does not affect the amplitude of the excur- 
sion. The same value of a and the same expression (0} 
for /(s) are expected for any distribution of single steps 
P(£) with finite variance, for which central limit theorem 
holds. Eq. (J2J can be easily computed explicitly in the 
case of a bimodal P(£). 

To extend our analysis to other processes we resort to 
numerical simulations for the evaluation of (x(£))t- We 
consider walks starting at x — at some negative time, 
and we take as t = 1 the first time such that x(t) > e, 
where e is a small positive threshold. We average over all 
walks that first return between [— e, e] at a specified time 




FIG. 2: Main: Plot of the normalized average fluctuation 
{x(t))r for a random walk with power-law distributed steps 
with jj, = 0.7. The solid line represents the universal func- 
tion @. Small circles are for several values of T ranging from 
199 to 9699. Inset: Plot of the normalization factor N(T), 
which scales as T 1 ^, indicating that a — 1/fj,. The dashed 
line is proportional to T 1 '^. 



T, under the constraint x(t) > e for 1 < t < T, For each 
time T the average shape is normalized by the factor 
N(T) = J* ds(x(sT)) T . If the scaling form © holds, 
the normalized shapes for different T collapse on a single 
curve and N(T) grows as T a . We average from 10 4 to 10 7 
pulses for each T. However, as shown in experiments 0], 
smaller statistics (~ 10 3 pulses) is already sufficient for 
reasonably clean curves. 

it) Levy flights: While it is well known that ran- 
dom walks with finite variance of the step distribution 
belong to the Gaussian universality class, distributions 
with fat power-law tails give rise to a completely dif- 
ferent behavior |12| . For this reason we consider now 
process Q where the distribution of independent incre- 
ments £ has infinite variance, i.e. with tails decaying as 
p {€) oc \i\~^~ x where < fi < 2. Although (x(t)) T 
can still be written in the form ©, F(x,t) can not be 
computed by the image method. We find numerically 
that also in this case the scaling form 1S obeyed, 
with a = 1/ ji and the scaling function /(s) of the same 
semicircular shape Q of the case with finite variance. 
This holds for any fi > (in Fig. |2]we show {x(t))x for 
fi = 0.7), thus even regardless the existence of the first 
(absolute) moment of the single step distribution 0] . 

While a scaling exponent a equal to the wandering ex- 
ponent of the free process could be expected [l4| , the fact 
that the shape of the scaling function remains exactly the 
same is surprising. Note that the short time expansion 
of Eq. © 



(x(t)) T oc T 1 ^- 1 / 2 t 1 ' 2 



(5) 



in the Gaussian case fi = 2 gives (x(£))t oc t 1 / 2 , equal 
to the wandering of the unconstrained walk. One could 
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FIG. 3: Average excursion for a random walk in a parabolic 
well [Eq. (J] for 1/A = 20. From top to bottom lines are for 
T = 10, 25, 50, 100, 250. Notice that the shape flattens to the 



constant value y/i/J^rX) (thin line). 



naively expect by analogy a short time behavior as t 1 ^ 
for fj, < 2 (corresponding to the free walk). Eq. (JSJ shows 
that this is not the case, indicating that the Levy flight 
feels the constraints even for short times. 

Universality in the shape of the scaling function is not 
restricted to changes of the variance (£ 2 ) of the single 
step distribution. Also variations of the mean value (£), 
corresponding to the addition of a bias to the walk, do 
not modify f(s). This is easily shown analytically in the 
case of gaussian or bimodal P(£) and has been checked 
numerically for the other cases. Therefore we can con- 
clude that, for all processes of the type (Q with uncorre- 
cted noise, the scaling function f(s) is proportional to a 
semicircle. 

Hi) Walk in a parabolic well: A further generalization 
is to consider a walk in a harmonic potential 



x(t + 1) = x(t) - Xx(t) + f (i) 



(6) 



where the term Ax describes the effect of a parabolic 
well. The damping A introduces a characteristic time, of 
order 1/A, so that scaling breaks down for large times. 
The probability for the walker to return at time T, i. e. 
the distribution of return times P(T), remains the same 
(oc T~ 3 / 2 ) for T <C 1/A, while it decays exponentially for 
T 3> 1/A. Also for this process the average fluctuation 
can be computed analytically via the image method and 
is 11] 



(x(t)h 



4 (1 - e- 2At )(l - e -2A(T-t)) 



ttA 



1 



,-2AT 



(J) 



As shown in Fig. |3 for T <C 1/A we get a semicircle, 
the same result as for A = 0. For larger T instead, as 
correlations decay, the form of the fluctuation tends to 
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FIG. 4: Main: Plot of the normalized average fluctuation 
{x(t))T for a random walk fed by noise with short time mem- 
ory (r = 1000). Empty symbols are for uniformly distributed 
n, while the filled ones are for power-law distributed n with 
fi = 3/2. Circles are for a time of the order of r, while tri- 
angles are for T 3> t. The solid line represent the universal 
function after normalization. Inset: The tails for small t/T 
of the same curves of the main panel plus a straight dashed 
line of slope 1. 



become flat, while keeping the (t/T) 1 / 2 ((1-i/T) 1 / 2 ) be- 
havior at the left (right) tail. This is strongly reminiscent 
of the flattening observed experimentally by Durin and 
Zapperi in the study of the Barkhausen noise in magnetic 
systems |l5j : indeed, in that case the damping term is 
provided the demagnetizing field, which controls the cut- 
off in the avalanche distribution ^(| . 

So far we have considered processes with increments 
extracted independently at each time step. Let us 
now analyze the effect of temporal correlations in the 
stochastic noise. More precisely, we consider processes 
such that 



g{t,t') = mm'))-mmt'))*Kt> 



(8) 



iv) Short-ranged memory: Let us first study the case of 
noise with correlations decaying exponentially over some 
interval r, i.e. g(t,t') — exp (— \t — t'\/r). The average 
fluctuation for this process, evaluated numerically, is re- 
ported in Fig. 01 showing the existence of two regimes. 
When T w r the shape of (x(£))t is parabolic, similarly 
to what has been observed in Q for the random field 
Ising model. The detailed form depends on the distribu- 
tion of single steps P(£) and in particular it is skewed 
for power-law distributed £. However, some form of uni- 
versality with respect to P(£) is present also in this case: 
The tails of (x(t))x, which go to zero linearly both for 
small s [Fig. 01 inset] and for s — > 1 (not shown). In the 
long time limit r>T, memory in the noise is lost and, as 
expected, we recover the semicircular shape (Fig.0J), that 
characterizes processes with uncorrelated increments. 
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v) Long-ranged memory: We finally analyze a stochas- 
tic process with long-ranged memory, given by Eq. IjTjI. 

where now = J2i=i t an< ^ ^(0 ^ s an uncorre- 
lated noise. This process can be seen as a temporally 
discrete version of the continuous process dfx = r)(t), 
which is known as the random accelerated particle (RAP) 
and has been studied in the context of inelastic col- 
lapse of granular matter ^t|. For a distribution of n(t) 
with unitary variance the correlation function of £(t) is 
g(t,t') = min(t,t') and hence it does not decay when 
the difference between t and t' grows. The first return 
probability to the origin of the walker has been recently 
shown analytically to decay as P(T) oc T~ 5 / 4 0. By 
evaluating numerically the average fluctuation (x(t))x we 
recover also in this case the scaling form [Eq. J3J], with 
a = 3/2, that again corresponds to the wandering ex- 
ponent for the unconstrained process. The normalized 
scaling function, plotted in Fig. is indistinguishable 
from the simple curve 



x 3 / 2 {l-x) 
5(3/2,1) 



(9) 



where B is the Beta function. For a RAP with power- 
law distributed rj (infinite variance), the exponent a is 
1 + l//i (again as in the free case). As in the case with 
short-ranged correlations the scaling function differs from 
the case with finite variance (Fig.[SJ, and the asymmetry 
is more evident for small /i. Nevertheless the behavior of 
the tails remains universal, being /(s) oc s 3 / 2 for s — ► 
(see Fig. inset) and f(s) oc 1 — s for s — > 1 (not 
shown) . 

In conclusion, we have tackled theoretically the prob- 



lem of analyzing generic time series from the point of 
view of the average shape of their fluctuations. In par- 
ticular, we have studied this quantity for stochastic pro- 
cesses of the form {Q. We find that the scaling law 
{x(t))r — T a f(t/T) is obeyed in an appropriate time 
regime. The exponent a coincides with the wandering ex- 
ponent of the unconstrained process, as can be expected 
on the basis of scaling arguments. More interestingly, 
the scaling function /(s) exhibits some degree of uni- 
versality. For any process with uncorrelated increments 
f(s) is proportional to a semicircle, with no dependence 
whatsoever on the distribution of single steps. Modifi- 
cation of this shape can be obtained only via additional 
x-dependent terms in Eq. (JTJ or introducing correlations 
in the stochastic increments. In this last case, the aver- 
age fluctuation is in general asymmetric and, although its 
detailed shape depends on the distribution of single in- 
crements, the tails are universal, depending only on the 
short- or long-ranged nature of noise correlations. Be- 
cause different universal behaviors encode solely correla- 
tion properties of the signal, this quantity is a powerful 
tool in the statistical analysis of temporal series: In many 
systems it would be very interesting to compare quanti- 
tatively (x(t))x between experiments and models in the 
light of our results. In the case of Barkhausen noise, 
where this comparison has pointed out the limitations 
of current models, our study provides a guide for refin- 
ing them in order to better capture the physics of the 
phenomenon. 
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FIG. 5: Main: Plot of the normalized average fluctuation 
(x(t))T for the random accelerated particle. Triangles are 
for uniformly distributed rj, while squares and circles are for 
power-law distributed rj, respectively with /i = 1.5 and /i 
0.7. The values of T range from 1000 to 100000. Inset: The 
tails for small t/T of the same curves of the main panel plus 
a straight line of slope 3/2. 
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